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We evaluate the thermal corrections to the generalized Gaussian effective 
^ potential. We carry out the calculations of the lowest order corrections in the case 

of self-interacting scalar helds in one and two spatial dimensions, and study the 

O 

l/-} restoration of the symmetry at high temperatures. 
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In a previous study [1,2] we introduced the generalized Gaussian effective 
potential which allows to evaluate in a systematic manner the corrections to the 
Gaussian approximation. In this paper we would like to study the thermal correc- 
tions to the generalized Gaussian effective potential. In particular we are interested 
in the phenomenon of the symmetry restoration at high temperatures. To this end 
we shall discuss in details the lowest order thermal corrections for self-interacting 
scalar fields in one and two spatial dimensions. 

Let us firstly discuss the thermal corrections in the general setting of self- 
interacting scalar fields in v + 1 dimensions. In the fixed-time Schrodinger repre- 
sentation, the Hamiltonian reads 



H= / d v x 



(1) 
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As it is well known [3], the Gaussian effective potential is given by 

VcepUo) = v rmn\<, 0> < ^ > , (2) 
where V"is the spatial volume and \tpo > is a Gaussian trial functional centered at 

The generalized Gaussian effective potential is set up as follows. Starting from 
the vacuum functional \tpo >, one constructs a variational Fock basis \n >. After 
that, the Hamiltonian H is splitted into the unperturbated Hamiltonian Ho and 
the perturbation Hj. Ho and Hi are defined respectively by the diagonal and 
off-diagonal elements of H on the basis \n >. By switching the interaction Hi on, 
the unperturbed ground state \tpo > evolves into \Q >. Whereupon the generalized 
Gaussian effective potential is given by [1]: 



(3) 



with the constraint 



-o 

dtr 



<"!'(*>!«> = ^ (4) 

Using the Gell-Mann and Low theorem on the ground state [4], we showed that 

i 00 (—{) n r° r 
V G (<f> ) = V GE p(4>o) + V Y1 / dtl - / 

v n =i n ' J -°° J ~ 

■ < 0\T(H I (0)H I (t 1 )...H I (t n ))\0 > con nected e^ + - +t ^, (5) 

where |0 >= \tpo >, and Hi(t) is the interaction Hamiltonian in the interaction 
representation. 

We, now, turn to finite-temperature field theory. Thermal equilibrium is achieved 
when the system has minimized its free energy F. Standard thermodinamics [5] 
gives 

F=-\lnZ, (6) 

P 

where Z is the partition function 

Z = tr (e-? H ) (7) 
and (3 — y (Boltzmann's constant k — 1). 

The finite temperature generalized Gaussian effective potential is defined by 

Va(<f>o) = y = -^lnZ (8) 

with the constraint (4). 

In order to evaluate the partition function we face with the problem of calcu- 
lating the trace in Eq. (7). Indeed, the trace involves a summation over the 



unknown eigenstates of the full Hamiltonian H . However, in our approach the full 
Hamiltonian has been splitted into the unperturbated Hamiltonian Ho and the 
perturbation Hj. Thus the natural strategy is to perform the trace in Eq. (7) 
by employing the thermodinamic perturbation theory [6]. Indeed, to first order 
in Hj, the trace can be replaced by a summation over the eigenstates of Ho. A 
straightforward calculation gives: 

Z = Z (l-p<H I > 13 ), (9) 

where 

Z = tr{e- fjHo ) = < n\e- fjHo \n > (10) 



and 



< Hi > p = ^trie-P^H!) = ^- ^ < n\e- pH °H r \n > . (11) 



n 



In Equation (10) and (11) we sum over the eigenstates of free scalar field Hamil- 
tonian with mass /i, where ji satisfies the gap equation: 



n n A , ? A f d"k 1 

MW + -^ + -/^^=. (12) 



Observing that 

<rc|i?/|rc>=0, (13) 
and that the states \n > are eigenstates of Ho, one readily gets 

< H T > fj = 0. (14) 
As a consequence, a straightforward calculation gives: 
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V£(*o)=V GB p(*o) + ^ J -^ln(l-e-^) ) + o[( Hi n (15) 



where g(k) — \jk 2 + /i 2 . Note that Eq. (15) differs from the finite-temperature 
Gaussian effective potential [7]. The difference resides in the different use of the 
gap equation. In our scheme the gap equation (12) is fixed once and for all. In 
particular it does not depend on the temperature . On the other hand, in the finite- 
temperature Gaussian effective potential approach the gap equation includes the 
thermal corrections. Now, the gap equation fixes the basis to sum over in the 
thermal average. This means that different gap equations lead to inequivalent 
basis. As a matter of fact the discrepancy between our result Eq. (15) and the 
finite-temperature Gaussian effective potential comes from the thermal average of 
the interaction Hamiltonian. In our approach Eq.(14) holds, whereas in Ref. [7] 
< H T 0. 

As we have discussed in Ref. [1,2] the main advantages of our approach at zero 
temperature reside in the diagrammatic expansion of the higher order corrections 
which is amenable to a diagrammatic resummation. Moreover the renormalization 
of the generalized Gaussian effective potential relies on the underlying </>o = field 
theory. Remarkably enough, it turns out that these fine features extend also 
to the case of scalar field theory at finite temperatures. Indeed we can write a 
manageable formula for the higher order corrections to Eq. (15). To do this, the 
usual thermodinamic perturbation theory is useless. Instead we may follow the 
Matsubara's methods [8]. In the Matsubara's scheme one deals with scalar fields 
which depend on the fictious imaginary time r varying in the interval (0,/3). 

If the Hamiltonian of the system in thermal equilibrium can be written as H — 
Ho + Hj, then one can show that the corrections to the thermodynamic potential 
are given by 
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Ail = -jjln < T T exp [ - \ H I (r)dT ) >, (16) 



where Hj(t) is the interaction Hamiltonian in the Matsubara's interaction repre- 
sentation, TV is the r-ordering operator, and the thermal averages are done with 
the free field partition function. Moreover, it turns out that only connected dia- 
grams contribute to AQ: 



I dTl - Tm < TriHiin)-^^)) >i nnected . (17) 



In our case, if we write 



V2(<h) = V GE p(<h) + ^ f ^n(l-e-^))+A^(^o), (18) 



we readily get 



1 00 (— \) m ff 3 

AV G (cf> ) = -- £ i-f- / dn...r m < T T {H I {n)...H I {r m )) >i nnected . (19) 



Note that, due to Eq. (14), the sum in Eq. (19) starts from m — 2. The 
thermal average of time-ordered products is evaluated by means of the Wick's 
theorem for thermal fields [9]. In this way we obtain the thermal corrections to 
the generalized Gaussian effective potential as connected vacuum diagrams. In 
Figure 1 we display the second order thermal corrections. The vertices can be 
extracted from the interaction Hamiltonian [1,2]: 



H T = / d v x 



(20) 
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where r/(x) = (f>(x) — <f>o- 

The solid lines in Fig.l are the thermal propagators of the scalar fields with mass 



Gp(x-y,n-T 2 ) =< T T ri(x,T 1 )ri(y,T 2 ) > p = - J 



I +^ f d^fc e i[k(S-y)-w n (T 1 -T 2 )] 

(21) 



where uj n — ^p. Note that the thermal propagator is periodic in the time variable 
with period 

A distinguishing feature of the graphical expansion of Eq.(19) with respect to 
Eq.(5) stems from the fact that the (ml) -1 coming from the rath order term is not 
completely cancelled by the number of different Wick contractions corresponding 
to a given graph. Consequently, a graph contributes to AVq (cf>o) in proportion 
to a combinatoric coefficient depending on the order of the graph. Moreover in 
evaluating the contribution due to a given graph one should take care of the normal 
ordering prescription in the interaction Hamiltonian. It turns out that the normal 
ordering in Hi modifies the so-called anomalous diagrams, i.e. the diagrams which 
vanish at zero temperature [10]. For istance, in Fig. 1 the diagrams (b) and (c) 
are anomalous. We find that the normal ordering of the interaction Hamiltonian 
modifies the diagrams like (b) and (c) in Fig. 1 by replacing G^O) with G ( g(0), 
where 



G /9 (0) = G /9 (0)- lim G /9 (0) = 

/3->- 00 



(2tt) 



2g(k) 



cotgh 



(5g(k) 



(22) 



A full account of the second order thermal corrections to the generalized Gaussian 
effective potential will be presented elsewhere [11]. In the remainder of the present 
paper we would like to analyze the lowest order thermal corrections, Eq.(15), for 
the case of v — 1,2 spatial dimensions. 
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In one spatial dimension Eq.(15) reads 



V£(* ) = V GE p{fa) + ff 2 J™ dt ln ( X " e"^+^) , (23) 

where f3 — /3/j,q,x — and /io = m(</>o = 0). In Figure 2 we show V^^o) — ^g(O) 
(in units of /ig) versus </>o for A > A c (*). As we can see, the symmetry broken 
at T=0 gets restored for T > T c . Obviously the critical temperature depends on 
A. It turns out that T c can be extimate, within a few percent, by means of the 
high-temperature expansion of the integral in Eq. (23). By following the classical 
work by Dolan and Jackiw [12] we find the following high-temperature expansion: 



- Vgep{4>q) + Mc 



— ^+ 

2/3 2 2/3 



jl + jL ln ( £2/1) _ m. 

87T 47T \ 47T / 647T 3 



x 2 f3 2 + 



C(5) 

5127T 5 



(24) 



where Q{z) is the Rlemann's zeta function. 

In Fig. 2 we also show the high-temperature expansion Eq. (24) (dashed line); 
we can see that the high-temperature expansion is a good approximation even at 
T ~ T c . Indeed For A = 4 we find T c ~ 1.27, which agrees with the value obtained 
by using the high-temperature expansion (24). 

The case of two spatial dimensions can be dealed with in a similar way. We have 



V G {® ) = V GE p{$o) + 



{2n)f3 3 Jo 



(25) 



,2 

where now <J>n — • 

The high temperature expansion of the integral in (25) gives [12] 



M Notations as in Ref. [2] 



8 



V£($ ) = V GE p(&o)+f4 



7T^ T 1 7*3/2 2 



. 90/3 4 24/3 2 12tt ^ 64tt 2 

ar 2 +0(/? 2 )l (26) 



64tt 



where the constant c = | + 2/n47r — 2y, and 7 is Euler-Mascheroni's constant. 
In Figure 3 we display Eq. (25) together with Eq. (26) in units of /ig (we subtract 
the temperature dependent constant Vj {§)). Again, the thermal corrections lead 
to the expected symmetry restoration at high temperatures. In this case, however, 
the high-temperature expansion Eq.(26) is a valid approximation for T » T c only. 
So that the critical temperature T c must be extimate by numerical methods. 

In summary, we have discussed the thermal corrections to the generalized 
Gaussian effective potential. In the lowest order approximation we obtain the 
restoration of the symmetry at high temperatures by a first order phase transition. 
The second order thermal corrections are under study. We plan to present the 
results in a future publication [11]. 
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FIGURE CAPTIONS 



Fig. 1 Second order thermal corrections to the generalized Gaussian effective poten- 
tial. 

Fig. 2 Thermal corrections to generalized Gaussian effective potential to the lowest 
order for v — 1 and A = 4. Dashed lines refer to high-temperature expansion. 

Fig. 3 Thermal corrections to generalized Gaussian effective potential to the lowest 
order for two spatial dimensions and A = 4. Dashed line refers to high- 
temperature expansion. The critical temperature is T c ~ 1.60. 
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